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Abstract 

The paper provides a rigorous analysis of the dispersion spectrum of SH (shear 
horizontal) elastic waves in periodically stratified solids. The problem consists of an 
ordinary differential wave equation with periodic coefficients, which involves two free 
parameters lo (the frequency) and k (the wavenumber in the direction orthogonal to 
the axis of periodicity). Solutions of this equation satisfy a quasi-periodic boimdary 
condition which yields the Floquet parameter K. The resulting dispersion surface 
uj{K, k) may be characterized through its cuts at constant values of K. k and uj that 
define the passband (real K) and stopband areas, the Floquet branches and the isofre- 
quency curves, respectively. The paper combines complementary approaches based 
on eigenvalue problems and on the monodromy matrix M. The pivotal object is the 
Lyapunov function A (w^, fc^) = ^traceM = cosK which is generalized as a function 
of two variables. Its analytical properties, asymptotics and bounds are examined and 
an explicit form of its derivatives obtained. Attention is given to the special case of 
a zero-width stopband. These ingredients are used to analyze the cuts of the surface 
uj{K, k). The derivatives of the functions w(/c) at fixed K and uj{K) at fixed k and of 
the function K{k) at fixed w are described in detail. The curves u!{k) at fixed K are 
shown to be monotonic for real K, while they may be looped for complex K (i.e. in 
the stopband areas). The convexity of the closed (first) real isofrequency curve K{k) 
is proved thus ruling out low-frequency caustics of group velocity. The results are rele- 
vant to the broad area of applicability of ordinary differential equation for scalar waves 
in ID phononic (solid or fluid) and photonic crystals. 
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1 Introduction 



The wave equation with periodic coefficients is ubiquitous in physics and engineering. Its 
apphcations in acoustics of solids have gained a new momentum since the introduction of 
artificial periodic materials such as phononic crystals. A common mathematical framework 
is the Floquet-Bloch theory of partial differential equations with periodic coefficients |16j . 
It does not however yield many explicit results for the general case of 2D or 3D periodicity 
and vector waves. The notable exception allowing an explicit analysis is the case of ID 
periodicity and scalar waves which is governed by Hill's equation [TTj. The spectral proper- 
ties of Hill's equation are very well understood for the situation where the wave propagates 
along some fixed direction (parallel to the periodicity axis or not). This case implies a 
single spectral parameter. The objective of the present paper is to take on a broader per- 
spective of arbitrary (2D) propagation of scalar waves in ID periodic media. This setup 
implicates dependence on two spectral parameters and thus leads to more elaborate wave 
spectral properties. The specific problem to be addressed is described next. 

Consider SH (shear horizontal) wave motion of the form Uz {x, y, t) = U{y) exp [i (kx — ujt)] 
which travels in the symmetry plane XY of a stratified monoclinic elastic solid with pe- 
riodic density p{y) = p{y + T) and stiffness Cij^i (y) = Cijki (y + T). The elastodynamic 
equation yields a second-order ordinary differential equation for the amplitude U{y), 

dj {cijkAuk) = pui =^ (C44C/' + ikci^U)' + ik (C45C/' + ikc^^U) = -pu'^U, (1) 

where di = d/dx, 82 = d/dy, ' = d/dy and Voigt's indices 4 = yz, 5 = xz are used [3]. 
It is convenient to pass from U to u = Ue^^ with Lp (y) = ik p (C45/C44) dy which reduces 
([T])2 to the Sturm-Liouville form 

{pi{y)u\y))' - k^^i2{y)u{y) = -u?p{y)u{y), (2) 

where pi = C44 and p2 = C55 — CI5/C44 denote the shear moduli. Equation ^ is the object 
of our study. The coefficients pifi{y) and p{y) are T-periodic strictly positive piecewise 
continuous functions of y € M, and k, u are two real parameters (unless otherwise specified). 
The functions u{y) and pi{y)u'(y) are assumed absolutely continuous. They satisfy the 
quasi-periodic boundary conditions 

u{T) = e^^^n(O), pi{T)u'{T) = e'^^ pi{0)u'{0) (3) 

with the Floquet parameter K ^ C, which by periodicity of e'^^'^ may be defined on the 
strip KeKT G [— tt, vr] called the Brillouin zone. Note that Eq. ([2]) admits equivalent 
representations obtained by changing the function and/or variable. For instance, replacing 
the variable y ^ y = /^^ /x^f ^ (?) d? recasts ^ in the form of a weighted Schrodinger 
equation 

u" (tj) + Lo^Z^u {y) = 0, with cj^Z^ = (w^ - /xs/cVp) ^O' ^0 = P/^i- (4) 
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Note that this transformation does not require reinforcing the above-imposed condition 
of piecewise continuity of /Ui(y). The coefficients Z and Zq {Z = Zq at k = 0) have the 
physical meaning of, respectively, impedance and normal impedance that we will find useful 
for interpretations. 

There exists a comprehensive spectral theory describing the eigenvalues io'^ [n € N) of 
([2]), (|3| as functions of K at fixed k, e.g. P ESI El [221 dSl [T] . From this perspective, the 
spectrum for real ii' G M is represented by the Floquet branches uJniK) on the (a;, K)-plane. 
Each branch spans a finite range on the cj-axis, called a passband, with a corresponding 
bounded solution Un{y)- Separating them are the ranges of w, called stopbands, where 
w G M and KT G vrZ + i (M\0) . Properties of the functional dependence of uJn{K) at 
fixed k can be described by various analytical means. One of the key ingredients of this 
theory is the so-called Lyapunov real- valued function A(ci;^) defined as the half trace of the 
monodromy matrix (the propagator over a period). By this definition, A(a;^) = cos KT 
determines the passbands and stopbands as the ranges |A(a;^)| < 1 and |A(u;^)| > 1, 
respectively. 

The present work is concerned with the more general framework in which the parameter 
k is considered as an independent variable on top of u and K. Keeping as an eigenvalue 
of Eqs. ([2])-([3]) now implies its dependence on two parameters: Un = ujn{K,k). For K 
real, uJn{K,k) is a multisheet surface whose sheets projected on the (w, /c)-plane span the 
passband areas bounded by the cutoff lines (| A| = 1) and separated by the stopband areas. 
Cutting this surface by the planes k = const and uo = const produces the Floquet branches 
and the isofrequency (a.k.a. slowness) curves, respectively. Clearly, such perspective is 
considerably richer than the one restricted to the Floquet curves at fixed k. It is also 
important to note that the present study differs from the two-parameter Sturm-Liouville 
problem with Dirichlet, Neumann and Robin boundary conditions, which has been studied 
elsewhere, see e.g. [UIMj- 

The structure and main results of the paper are as follows. Section 2 introduces com- 
plementary approaches based on differential operators Axik), Bk{^) defined by ([2]), ([s]) 
and on the matricant M {y, yo) of the equivalent differential system. The operators Axik), 
Bk{^) are self-adjoint and have a complete orthogonal system of joint eigenfunctions, as 
shown in Appendix Al by explicit construction of their resolvent operators. The eigenval- 
ues and k"^ of Axik) and I3k{^) are then linked to the monodromy matrix M(T, 0) 
with eigenvalues e^^^ via the generalized (depending on two parameters) Lyapunov func- 
tion A{uj'^,k'^) = ^traceM(r, 0) = cos KT. Section 3 describes this function in some detail. 
It is shown in §3.1 that A(t<;^,A:^) inside the passbands |A| < 1 has non-zero first deriva- 
tives in both and k'^, and that A(a;^) for fixed k'^ and A(fc^) at fixed each satisfies 
Laguerre's theorem (by virtue of the estimates of A(w^, k"^) given in Appendix A2). These 
two fundamental facts explain the regular structure of the passband/stopband spectrum 
on the (w, A;)-plane. The WKB approach [lOj is used in §3.2 to provide an insight into the 
asymptotic behaviour of stopbands for continuous and piecewise continuous periodic coef- 
ficients. Zero-width stopbands (ZWS) are introduced and analyzed in §3.3. Generalizing 
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the concept of degenerate gaps of a one-parameter spectrum (e.g. [191 [131 E]), ZWS are 
intersections of the analytical cutoff curves |A| = 1 with the (a;, A;)-plane. It is shown that 
ZWS may or may not exist for an arbitrary periodic profile of p{y) and /ii,2(y)) are likely 
to exist for any profile that is even about the period midpoint, and always exist for a peri- 
odically bilayered structure. In the model cases, ZWS may also form infinite lines on the 
(cj, /i;)-plane. Closed-form expressions for the partial derivatives of A(a;^, k"^) are obtained 
in §3.4. The derivative of any order is a multiple integral of the product of, specifically, right 
off-diagonal elements M2 of the matricant M taken at different points y within the period 
and weighted by p{y) and/or ^2(y)- An alternative representation is derived for the first- 
order derivatives of A(a;^, k"^) within the passbands by using the eigenfunctions of Axik) 
and Bk{^^)- The two equivalent formulas obtained for the first derivatives of A(a;^,/c^) 
provide an explicit meaning to their sign-definiteness and offer useful complementary in- 
sight. In particular, it reveals some interesting attributes of the function M2{y + 1,?/), 
whose zeros (w. A;) are y-dependent solutions of the Dirichlet problem on [y, y + T], see 
§3.5. The properties of the Lyapunov function A(a;^,fc^) (= cos KT) and the expressions 
for its derivatives established in Section 3 are then used in Section 4 to analyze principal 
cuts of the dispersion surface ujniK,k). In §4.1, dependence uj{k) for fixed K is studied. 
It is shown that if K is real then the curves 0Jn{k) are monotonic (this may not be so for 
complex K) and they tend to the same linear asymptote A; min^^gjo.T] [/^2(y)//'(y)] which 
is independent of n. In §4.2, the dependence io{K) at fixed k is discussed. For real K, 
the first non-zero derivative of Floquet branches ujn{K) is provided (it is a first derivative 
inside the passbands and a second one at the cutoffs); for the stopbands, the condition on oo 
realizing maximum of |Imi(' (ti;)| is formulated. The real isofrequency curves K{k) at fixed 
uj are considered in §§4.3 and 4.4. Particular attention is given to the closed isofrequency 
curve arising for u) less than the first cutoff wi (7rT~^,0) . It is proved that, whatever the 
distortion of its shape due to unidirectional periodicity may be, this isofrequency curve is 
always convex and hence low- frequency caustics of the group velocity Vw are impossible. 
Finally, useful bounds on the first eigenvalue LOi{K,h) for KT S [— 7r,7r] and any k are 
provided in Appendix A3. 

Without loss of generality, in the following we take T = 1; more precisely, this implies 
the redefinitions y =^ y/T = y, oj ^ ojT = u, k ^ kT = k and K ^ KT = K so that the 
variables y and w, k, K are hereafter non-dimensional. We also assume throughout that 
T = 1 is a minimal possible period. 

2 Eigenvalue problem, monodromy matrix and Lyapunov 
function 

Equation ^ with the conditions ^ can be considered in either of the equivalent forms 

Aru = w^n, Bku = k^u, u G Dk (5) 
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with the operators Ak ^ Axik) and Bk = Bk{oj) 

Aku = — (uiu')' + k'^—u, Bku= — (uiu)' + u}'^—u. (6) 
P P M2 ' P2 

Their common domain is 

DK = {ueD: 7^il) = e^^Vm}, ( u{y) \ 

D = {u(^AC%i\, /.in'G^c[o,i]}, "^^y^ Vw(yK(y)y' 

where K G C and AC[0, 1] is the space of all absolutely continuous functions from [0, 1] to C 



(note that using "z" in the definition of rj and hence in ( 10 )2 is a conventional option which 
is useful for a compact form of (13)i and similar identities). Let {■■,-)p and ||-||^ be 



a standard inner product and norm in the Hilbert space Hp^ = (0; 1) of functions 

with quadratically summable measure p{y)dy and /i2(y)dy, respectively; so that 



{u,v)o= I P{y)u{y)v*{y)dy, \\u\\l = {u,u), 







1 

iu,v)p^ = 1^ fj.2iy)u{y)v*{y)dy, \\u\\l^ = {u,u)p^ 



(8) 



where * means complex conjugation. 

The operator ^ on (M) with eigenvalues co"^ (or k'^) can be represented as a direct 
integral decomposition ®K&[o.2n]'^K (or (Bk(^[o,2tt]^3k) [22j. Therefore the spectrum of the 
operator ([2]) is a union of all eigenvalues of Ak (or Bk) for K £ [0,27r] and hence for all 
K £ R. since Ak = Ak+2-k-, Bk = Bk+2'w The operators Ak and Bk are symmetric if 
G M, i.e. {Aku, v)p = {u, Akv)^ , {Bku, w)^^ — ^kv)^^ for n, f G Dk, and they both 
have compact and self-adjoint resolvents that satisfy the Hilbert-Schmidt theorem (see 
Appendix Al). Therefore Ak and Bk are self-adjoint with purely discrete spectra a (Ak) 
and a (Bk) containing an infinite number of real eigenvalues a;^ {K, k) and k"^ {K, uj) (n € 
N), and corresponding eigenfunctions Un{= Un,A and m^.b) forming a complete orthogonal 
system in the spaces Tip and V.^^ , respectively. The operator Ak is positive for any /c G M 
(i.e. for any A;^ > 0), 

{Aku, u)p>0 (> at /c / 0) , (9) 

so its spectrum a (Ak) consists of non-negative eigenvalues LJ^{K,k) (strictly positive at 
k ^ 0), which are hereafter numbered in increasing order cji < UJ2 < ... By contrast, 
Bk is not sign-definite and hence its spectrum a (Bk) includes both positive and negative 
eigenvalues k'^{K,u:). Note that real eigenvalues of Ak and Bk are also admitted at 
luiK 7^ (see Definition Qc) below). 
Equation ^ can be recast as 

r/'(2/) = Q(y)r?(2/) withQ(y)=if ^2°_^^2 " ') (10) 
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ioi r}{y) introduced in Given an initial condition r/ (j/o), Eq. (10)i has a unique solution 



riiy) = M{y,yo)i]iyo) 
defined through the propagator matrix, or matricant, 



(11) 



M(y,yo) 



Mi{y,yo) M2iy,yo) 
M3(y,yo) M4(y 



^yO)J J yo 

I + r Q (<?i) d?i + r Q (?i) d?i Q (?2) d?2 + . . . 

J ?;n J «() J Jin 



(12) 



where / is the multiplicative integral evaluated by the Peano series [2lj and I is the 2x2 
identity matrix. Note that det M (y, yo) = 1 due to trQ = 0, where tr means the trace. By 
(|l0|) Q = -TQ+T for uj'^, k"^ and so 

M-i(y,yo) = TM+(y,yo)T ImMi,4 (y, yo) = 0, Re M2,3 (y, yo) = 0, (13) 

where ^ denotes Hermitian transpose and T is the 2x2 matrix with zero diagonal and unit 
off-diagonal elements. If Q (y) is also even about the midpoint of the interval [yo, y] then 



M(y,yo) = TM^(y,yo)T ^ Afi (y, yo) = M4 (y, yo) 



(14) 



where denotes transpose. The properties (13)i and (14)i are actually valid for matrices 
Q and M of arbitrary n x n size (see 



for details), while (13j) and (14)2 are attributes 
of the 2x2 case which admits easy direct proofs (e.g. (|13[)2 is evident from the definition 
([7])2 of r/ with a real scalar u). 

Assume a periodic Q(y) so that Q(y) = Q (y + 1) and hence M (y, yo) = M (y + 1, yo + 1). 
The propagator M (yo + l,yo) over a period [yo,yo + 1] is called the monodromy matrix. 
For any yo = y, denote its elements as 



M(y + l,y) 



mi(y) im2{y) 



mi,4(y) = Mi,4(y + l,y), 

im2,3{y) = M2,3{y + i,y), 



(15) 



where Im mj(y) = 0, j = 1..4, for w^, A;^ G M by (13)2. The assumed periodicity with use 
of the chain rule implies the identity 

M(y + l,y) = M(y + l,l)M(l,0)M(0,y) = M(y,0)M(l,0)M-i(y,0). (16) 

Remark 1 The trace and eigenvalues o/M(y + l,y) are independent of y by virtue of 

Hereafter, unless otherwise specified, we set yo = and define the monodromy matrix as 
M(1,0) with respect to the period [0, 1] (as in 
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Bearing in mind det M = 1, denote the eigenvalues of M(l, 0) by q and q ^. Introduce 
the generalized Lyapunov function 
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Aiuj^k^) ^ -trM{l,0) = -{q + q-') 



(17) 



which is analytic in uj by pOh, P2P and real for uj% /c^ G R by e3h- As noted 



above, the function A(a; , A: ) is independent of the interval on which the unit period is 
defined. It is also invariant for any similarity equivalent formulation of the system matrix 
Q(y) = C~^Q(y)C because trM = tr (C-^MC) = trM, leaving A(a;2,/c2) unchanged. 

Proposition 2 For any complex numbers k, uj, K, the following statements are equiva- 
lent: (i) uj'^ is an eigenvalue of the operator Axik); (ii) k'^ is an eigenvalue of the operator 
Bk{'^)', (Hi) k, Lo and K are connected by the equality 



A(w^A;2) 



cos K = 0. 



(18) 



Proof. The link (i)^(ii) follows from Eq. Consider (i),(ii)^ (Hi). According to (i) 
or (ii), up' or /c^ is an eigenvalue of, respectively, Ak{}^) or Bk{^)- Then there exists 
u{y) G Dk that satisfies ^ hence Q, and consequently the vector ??(y), generated by 
u{y) according to is a solution of Eq. So, by ([TT]), r/(l) = M(l,0)r/(0). On the 

other hand, as indicated in u{y) G Dk implies that r/ (1) = e*^r/ (0). Hence e*^ is an 
eigenvalue q of M(1,0), and the function A defined by (|17[ ) satisfies (18), that is (Hi). Now 
consider (iii)^(i),(ii). From (18) and the definition (17), the eigenvalue q of M(1,0) is 
q = e*^, and corresponding eigenvector w exists such that M(l,0)w = e*^w. Let u{y) be 



the first component of the solution r/(y) = M(y, 0)w of Eq. (10) with the initial condition 
77(0) = w. From the above, u{y) belongs to Dk and satisfies Eq. ([5]), which implies (i),(ii). 



Corollary 3 Each eigenfunction u of Ak and Bk is equal to the first component of the 
vector r](y) = M(y,0)w, where w is the eigenvector o/M(l,0) corresponding to the eigen- 
value q 



AK 



Definition 4 Passband areas, cutoffs and stopband areas are defined for uj'^ , A;^ € M (and 
hence real A{uu'^, k"^)) as follows: 

'|A|<1 {-i?^ K ^M.) passbands, 
{u, k) : <A = ±1 {^K e irZ) cutoffs, 

J A| > 1 {-^ K eirZ-^i (M\0)) stopbands. 

Before discussing general properties of the Lyapunov function A(a;^, k'^), it is expedient 
to mention its explicit properties at cj = and/or A; = 0. Obviously dA/duo = at w = 
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and aA/9A; = at A; = 0. By ^ and (fTTI 



A{u;^ k') = l + l (/.r'> ((/i2) - {p) uj') + 0((a;2 + k^) with (•) = / (•) dy; 

Jo 

dA/d{uj^) = -\ (p) (/xr^> , 5A/a (A:2) = \ (^-1) (^2) at = 0, = 0, 



(19) 



where the identity d? /g^' [/i (?) /2 (?i) + h (?) /i (?i)] d?i = (/i) (/2) was used in (19)i. 
Note that A(0, k^) > 1 for A;^ > and \dAld{u?)\^^^ < for k^ > 0, whereas the bounds 
of A(a;2,0) and the sign of [dA/d {k'^)]i^_Q are not fixed for > 0. Also note the explicit 
non-semisimple form of the matrix 



M(2/,0) 



1 -z/^vr'Wd? 



at cj = 0, /c = 0. 



(20) 



3 Properties of the Lyapunov function A(aj^,/c^) 
3.1 Formation of the passband/stopband spectrum 

We proceed with some observations on the analytical properties of the function A(a;^, k'^) 
that underlie the alternating structure of the passbands and stopbands. 

Lemma 5 Ifu^Rork^^R then A ^ [-1, 1] . 

Proof. If A S [—1,1] then according to Proposition [2] the identity (18) holds for G R 
and hence or /c^ is an eigenvalue of Axik) or Bk{(^), respectively. It was shown (see 



(11) and below) that the eigenvalues of Axik) are positive and the eigenvalues of Bk{^) 



are real. ■ 

Proposition 6 The derivatives dA/d{Lo'^) and dA/d(k^) do not vanish within an open 
passband interval A(a;^, k'^) G ( — 1, 1) . 

Proof. By Lemma [sj if A G (-1,1) then u}^,k^ G M. Suppose that dA/d{uj^) = for 
some real value uj'^. Then, because A(a;^) (= A(a;^,A;^) at fixed k) is an analytic function, 
there exists complex in the vicinity of cj^ for which A(a;^) G (—1,1) . This contradicts 
Lemma pi and hence dA/ d (uP'^ 7^ 0. The same reasoning proves that dA/d [k'^) 7^ 0. 



Consequently, Eq. (18) at fixed uj^ > (or fixed real k'^) has only real and simple roots /c^ 
(or a;2) if cosK £ (-1,1). ■ 

Proposition [6] plays a pivotal role in explaining the origin of the Floquet stopbands by 
the following simple reasoning. Consider p{y), Pi,2{y) resulting from an arbitrary peri- 
odic perturbation of some reference constant values pQ and poi,02, so that AiuJ^^k"^) is a 
perturbation of Ao(i^^,A;^) = cosK with = j^^'^ — ^jj^^"^- Since the first derivatives 
of A(cj^,A;^) do not vanish within (—1,1) , the perturbed extreme values Aq = ±1 must 
either remain equal to ±1 or exceed the range [—1, 1] , thereby leading to complex values 
K G ttZ + i (M\0), i.e., to the stopbands. 
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Proposition 7 For cj^, k"^ G M, the derivatives of any order n £ N of the functions A(cl!^) 
and A{k'^) (= /S.{uP' , k"^) at fixed k and fixed uj, respectively) have only real and simple zeros, 
each lying between consecutive zeros of the (n — l)th derivative of the same function. In 
particular, the first derivatives of A(a;^) and A(A;^) have a single and simple zero between 
consecutive zeros of A{uj'^ , k'^) and do not vanish elsewhere. 

Proof. It is shown in Appendix A2 that the functions A(cij^) and A(/c^) are entire functions 
of order of growth ^. Their zeros are the eigenvalues of the operators A.,^/2{k) and Bt^/2{^)^ 
and are therefore real and simple. Hence both functions satisfy the conditions of Laguerre's 
theorem (e.g. |27j). implying that the derivatives of A(a;^) and of A(A;^) are also entire 
functions with order of growth | and they have the desired properties. ■ 

Propositions [g] and [7] define the basic form of the function A(a;^, /c^) at fixed to or k. It 
is exemplified in Fig. [T] for a piecewise continuous profile of material coefficients chosen as 

/ii(y) = ^2(y) = ^(1 + 3y)'(2 + y), p{y) = 2 + y for y E [0, 1] (21) 

(taking /xi^2 in GPa and p in g/cm^ implies ujT = lo in MHz-mm in this and subsequent 
figures). Note that A((x>^) has an infinite number of zeros that are strictly positive and 
move rightwards as k increases, whereas A{k'^) has an infinite number of negative zeros at 
LO = which move one by one on the positive semi- axis A;^ > as a; increases. 




Since zeros of the first derivatives of A(a;^,A;^) cannot be points of infiection or zero- 
curvature by Proposition^ we can now refine the numbering of branches ujn{K,k) = 
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^ij:\{K^ k) (> 0) in the passbands as fohows: 

< oji{K, k) <uj2{K,k) < ... if K eR, K ^ nZ; 

<uji{0,k) <LO2{0,k) <LO3{0,k) <LOi{0,k) < ... if K G 27rZ; (22) 

< u;i(7r, k) < uj2{Tr, k) < uj3{tt, k) < tdiiir, k) < . . . if K £ n + 2txTL. 

With reference to (19) and Proposition [6| the sign of first derivatives of A(aj'^,A;^) along 



u:n{K,k) in the nth open passband |A| < 1 (see (([22|i)) is 

sgn [0A/9(w2)j ^ _ggj^ \^Q^iQ (^2-)j ^ _ ^23) 

The possibility of equality of two cutoffs (see (p2|)2,3), i-e. of a double root of the equation 
A(cl!^) = ±1, implies a zero-width stopband addressed in detail in §3.3. 

For the future use, let us also mention some properties of the Dirichlet and Neumann 
eigenvalues „ and „ of (p| satisfying the conditions u (0) = 0, it (1) = and u' (0) = 0, 
u' (1) = 0, respectively. It is Known that (jJy) n 

and a;N,n are simple zeros of the functions 
M2(l,0) and M3(l,0) of uj, which occur once per each stopband complemented by cut- 
offs (except the first stopband devoid of a;D,n)- The branches a;D,i(A:) < ojd^2 (k) ... and 



wn,i(^) < wn,2(^)--- are thus related to the passband eigenvalues ujn {K, k) of (22) as 



ujB,2j{k),u)N,2j+iik) G [u)2j{0, k),u)2j+i{0, k)]; u:B,2j-iik),ujt^,2j{k) G [a;2j-i(7r, A;), a;2j (vr, A;)], 

(24) 

where j G N and uj-^^i{k) G [0, wi(0, k)]. Recall that the stopbands and cutoffs are invariant 
with respect to the choice of the period interval [yo, yo + 1] = [0, 1] (see Remark[l]); however, 
the branches (^D,n(A;) and (^N,n(A;) within this area certainly depend on the choice of the 
point yo = 0. In other words, some fixed values uj, k realize the Dirichlet or Neumann 
conditions at the edges of [yo, yo + 1] iff yo is a zero of the function M2 {y + l,y) = im2{y) 
or M3(y + l,y) = im3{y), respectively (see §3.5 for further discussion). According to 



(14), if Q(y) is an even function about the midpoint of the period [yo, yo + 1] for some yo, 
then the Dirichlet and Neumann branches UB^nik) and uji^^nik) satisfying m2(yo) = and 
"^3(2/0) = coincide with the cutoff curves. We note the useful identity m2{y)m3{y) > 
for I A| < 1 which may be proved as follows: it obviously holds for A = due to det M = 1, 
and hence for any |A| < 1 due to the fact that m.2(y) and m2,{y) are strictly non-zero inside 



the passbands by (24). 



3.2 WKB asymptotics of A 

Some insight into the high-frequency spectrum in the case of continuous and piecewise 
continuous periodicity can be gained from the WKB asymptotics [10] of the Lyapunov 
function A(a;^,/c^) at fixed k. To this end recall the impedance Z = Zo\/l — \X2k'^ I P'^'^ 
with Zo = -s/pfJ-i introduced in Q. For any fixed k, let > k"^ maxy^^o,!] {P'2/ p) so that 
Z{y) is real (the so-called supersonic regime). Suppose for brevity that the overall periodic 
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profile of Z{y) has at most one point of discontinuity per period. If so, the zero-order WKB 
approximation A^^Lg of A takes an especially simple form 



A. 



(0) 

WKB 



[Z]l/2 + [Z]-^/")cOS {UJ 



1/2 ■ 



'Zdy), 



(25) 



where ±iu}fi^^Z are the eigenvalues of the matrix Q defined in ( 10 )2 and [Z] = Z [y^) /Z (yj") 
with Z {y^) = lim£_>o Z {yd ± e) is the relative jump of Z at the possible point yd of its pe- 
riodic discontinuity. Assume first that Z{y) is strictly continuous for any y (not restricted 



(0) 



WKB 



< 1 and thus can estimate zeros 



to [0, 1]) and hence [Z] = 1. Then Eq. (25) yields A 

of A but not the stopbands | A| > 1, whose widths (the frequency gaps between cutoffs, see 
( 19 12,3) may well be nonzero at finite co. Thus if Z{y) is continuous then Eq. ( 25 ) merely im- 



plies that the stopband widths tend to zero at any fixed k as u tends to infinity. The latter 
conclusion is also valid even if H2/P has periodic jumps but pfj^i is continuous throughout, 
so that [Z] 7^ 1 indicates existence of nonzero stopbands at finite uj but [Z] — t- [Zq] = 1 at 



a; — )• 00. On the other hand, if pfii does have a jump and so [Zq] / 1, then Eq. (25) shows 
that the stopband widths remain nonzero as cj — )• 00. Having stated this, we hasten to add 
that a physically sensible profile model should be related to the frequency u in that a finite 
w implies that a probing wave " sees" appropriately abrupt variations of material properties 
as jumps, which are of course smoothed out by the 'infinite zoom' of the limit u — )• 00. The 
above WKB conclusions on the high-frequency trends of cutoffs agree with a less general 
framework of, specifically, small periodic perturbations that provides expressions for the 
stopband widths through the Fourier series coefficients, see [31 [6]. 

As an example, consider again Fig. [l| which is plotted for a piecewise continuous 
profile (21) that gives [Z] = Vl-sJ (1 — Ah? /oj"^) / (4 — /uj"^) (note that a 'single periodic 
discontinuity yd is located at the edges of the period T = 1 by ( pTj ); however, similarly 
to Remark 1, A^^j^g does not depend on the choice of the period [0, 1] relative to yd)- It 
is easy to check that the exact curves A shown in Fig. [T^ are well fitted by the WKB 
approximation (25) (not displayed to avoid overloading the plot) once uj is greater enough 
than k max ^ ^2! P = 2A;. It is also seen from Fig. [l^ that increasing u) makes the curves 
A for different fixed k tend to that related to A: = 0, as predicted by Eq. (25). 



In the case of two or more discontinuity points per period, applying the WKB asymp- 



totics separately along each range of continuity modifies ( 25 ) to the form with two or more 
phase terms corresponding to the reflection-transmission at each discontinuity. For more 
examples of using the WKB approach to the periodic profile, see [23] . 
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3.3 Zero-width stopband 

3.3.1 Complementary definitions of ZWS 

The following definition of a zero-width stopband (ZWS|^ is motivated by the possible 



occurrence of the second and third cases in ( 22 ) . 



Definition 8 If uj = oj2ni0,k) = a;2n+i(0,fc) or u = uj2n-i{T^, k) = U2ni'^,k) for some 
uj, k GM. and n G N, then this cutoff point {u, k) is called a ZWS. 

It is essential that the cutoff curves are analytic (as any uJn{K, k) with fixed G M is, 
see §4.1), hence if two of them meet at a point they cannot conjoin. Thus an isolated 
ZWS implies intersection of two cutoff curves on the (w, A;)-plane and hence a saddle point 
|A| = 1 on the Lyapunov-function surface A{oj'^, k"^). For the same reason, if, exceptionally 
(see §3.3.3), a ZWS forms a line u) (k) of local extremum |A| = 1 of A(a;^, fc^), then such 
line cannot have an edge point. 

A comprehensive account of the properties of ZWS is based on the next proposition. 

Proposition 9 The following statements are equivalent: (i) {uj, k) is a ZWS; (ii) A(a;^, k"^) - 
±1 and dA{uj'^,k'^)/d{u^) = 0; (Hi) A{uj^,k'^) = ±1 and dA{uj'^,e)/d{k'^) = 0; (iv) 
M(1,0) = ±1. 

Proof. The link (i)<^(ii) follows from Definition [s] and Proposition[7j The link (i)^(iv) can 



be inferred e.g. via (24), which tells us that assuming (i) entails M2(l,0) = M3(l,0) = 



and hence Mi(l, 0)M4 (1, 0) = detM = 1, where Mi, M4 are real by ([T3j)2. Since (i) 
also means trM(l,0) = ±2, it follows that M(1,0) = ±1 as stated. Next let us show 
(iv)^(ii). Assume M(1,0) = ±1 for some w. A; e M. Note that A{uj'^ ,k'^) = ±1 by ( |l7| ). 
The (double) eigenvalue q = e*^ = ±1 of M(1,0) = ±1 has geometrical multiplicity 
2, hence uj'^ is an eigenvalue of Axik) of multiplicity 2 by Corollary [sj Now consider 
some K' G arbitrary close to K that yields cosK' = A{uP',k'^) G (—1,1). Since is 
a double eigenvalue of AK{k), the self-adjoint operator Ak'{^ has two distinct simple 

eigenvalues ^-fC', k^ close to u)'^ , and, by Propositions jij and jo these are distinct simple 

zeros of A(a;^,fc^) — cos A''. Therefore A(a;^,A;^) = ±1 is a local extremum of A(a;^,/c^), 
i.e. dA/d{ijj^) = at S^, A;^, which is equivalent to (ii). Note that reversing the above 



reasoning proves (ii)^(iv) without appeal to (24), and that invoking Bk{^) in place of 



Axik) provides a similar proof of (iii)<^(iv) (see also Proposition 16 below). 



Note that the point a; = 0, fc = which yields A = 1 is not a ZWS since it does not satisfy 



any of the above statements, which is evident from (19)-(20). 

Proposition [9] implies that the multiplicity of w^, k)^ as the roots of equation A(a;^, k"^) — 
cosK at K G M is the same as their multiplicity as the eigenvalues of Axik), I3k{^) (this 



^It is understood that a ZWS is actually not a 'stopband' (in the sense of Definition [4]). Note that a 
similar notion of 'zero-width passband' is inconceivable due to Proposition [7] 
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multiplicity is 2 at a ZWS and 1 elsewhere). This is noteworthy since such a parity does 
not always hold inside a 'true' stopband K ^M., where a double root or A;2 of Eq. (pl IS 



not a double eigenvalue of, respectively, Axik) or Bk{^) which are no longer self-adjoint 
for ^ M. It is also pointed out that the eigenvalue q = e*^ of M(1,0) has an algebraic 
multiplicity 2 at any cutoff, while its geometrical multiplicity is 2 only at cutoffs that are 
ZWS. 

Corollary 10 The matrix M(1,0) is non-semisimple for any cutoff {uj,k) unless it is a 
ZWS. 

We note that the non-semisimple nature of the monodromy matrix at the cutoffs has 
important ramifications for the interpretation of its matrix logarithm, which has been 
proposed as the basis for dynamic effective medium models, see pS} [26]. 

3.3.2 Considerations of the existence of ZWS 

To begin with, it is recalled that the period T = 1 is everywhere understood as a minimal 
possible period, so that trivial ZWS which turn up when T is a multiple of the minimal 
period are disregarded. 

Given an arbitrary periodic Q(y), the condition M(1,0) = ±1 stipulating existence 
of ZWS imposes three real constraints on two parameters to, k and hence is unlikely to 
hold. However, if the profile Q(y) is symmetric (even) about the midpoint of the period 



[0, 1], then, by virtue of (14), the above condition on M(1,0) implies only two constraints 
and thus such profile can be expected to yield a set of ZWS points (intersections of cutoff 
curves |A| = 1) on the {u,k)-plane. More precisely, since the cutoffs are independent of 
how the period interval is fixed (see Remark[T]), ZWS are expected to exist if a given profile 
Q{y) admits such a choice of the period interval [yo,yo + 1] = [0, 1] within which Q(y) is 
symmetric. 

Note that by definition any ZWS is also an intersection of Dirichlet and Neumann 



branches (24) while the inverse is generally not true. Moreover, in contrast to ZWS, the 
Dirichlet and Neumann branches and hence their intersections \^uj,k^^_-^ depend on the 
choice of the period interval. For instance, let Q(y) be symmetric with respect to a fixed 
period [0, 1] . Then the Dirichlet and Neumann branches coincide with the cutoff curves 
and hence any intersection {^,k}j-^^^ is a ZWS (see e.g. Fig. 1 of [24]). However, if for a 
given Q(y) = Q(?/ + 1) the period is shifted so that Q{y) is not even about its midpoint, 
then a new set {w, includes but generally does not coincide with the (unchanged) 
set of ZWS. 

As a simple explicit example, consider a periodically bilayered structure where Q,{y) 
takes two alternating constant values within two layers j = 1,2 that constitute a period 
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[0, 1]. The monodromy matrix is given by the standard expression 



Mil 0) - i '^2 cos V^i - 1^ sin ■02 sin V'l ^ cos ^^2 sin V'l + ^ sin tjj2 cos V'l \ 

' I iZ2 cos V'l sin ■02 + i-^i sin 0i cos 1^2 cos 02 cos "01 ~ §f sin 0^2 sin 0i j 

(26) 

where is the layer impedance defined in Q and tpj = ioZjdj/ fiij with dj for the 
layer thickness. The set of Dirichlet/Neumann intersections |a;,A;}j^_j^ is defined by 
simultaneous vanishing of both off-diagonal components of ( |26[ ), which implies the fol- 
lowing three options: (i) {sin^0i = 0, sin 0^2 = 0} , (ii) {cos0^i = 0, cos^02 = 0} and (iii) 
{Zi = Z2, sin (•01 + 02) = 0}, where (iii) may or may not hold for real lo, k [2]. It is seen 
that (i) and (iii) yield M(1,0) = ±1. Thus (i) and maybe (iii) define ZWS, while (ii) does 
not. 

Recall that an infinite periodically bilayered structure can always be considered over 
a three- layered period where the same stepwise profile Q(y) is symmetric. Hence the fact 
that any bilayered profile always admits ZWS (see e.g. Fig. 2b in §4.1) is consistent with 
the above conclusion that ZWS should be expected for the profiles Q(y) that can be defined 
as symmetric over some interval [yo ,yo + 1] . 

3.3.3 Model examples of regular loci of ZWS 

• Uniform normal impedance: Zq = p(y)/^i(y) = const at any y G [0, 1] . 

Let A; = 0. The coefficient in Q at A; = is Z (y) = Zq (y) , which is constant at Zo{y) 



const by virtue of fii > 0. Alternatively, note from (10)2 that Q(y) with k = and 
Zq = const has constant eigenvectors. Either of these observations readily shows that, for 
= 0, a dependence of w on > (not restricted to ii' G [0, vr]) is a straight line and thus 
all stopbands are ZWS, that is, there is no stopbands at all. The only difference with the 
case of constant p and pi is the slope of oj {K, 0) which is specified as follows: 

u{K,0) = KZo/{p) =K/Zo{iJ^^), (27) 

• Uniform speed: = P2{y)/ p[y) = const at any y G [0, 1] {pi{y) is arbitrary). 



The Lyapunov function is then A(cj , k ) = A(a; — c k , 0), from (10)2, and consequently 



Wn(i^, k) = yjujl {K,Q) + C^k^. (28) 
Hence if (-/rm, 0) with m = or 1 is a zero-width stopband, that is, if w„ (■/rm, 0) 



^n+i (Trm, 0), then by (28) a;„ (-/rm, k) = OJn+i (vrm, k) \/k, i.e. the entire line [uj"^ (vrm, k) , k) 



for any A; G M is a locus of ZWS. Note from (28) and (20) that the first cutoff (which is 
not a ZWS) is wi (0, A;) = ck = uj-^^i{k), where cjn,i(^) is the first Neumann solution for 

ye [0,1]. 
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Uniform normal impedance and speed: Zq = const and = const at any y £ [0, 1] . 



Now Eqs. (27) and (28) together imply that all stopbands are ZWS for any A: E M. Note 



that the inverse statement is true under an additional condition of absolute continuity of 
Zq, by the Borg theorem [5]. 

3.4 Explicit expressions for the derivatives of A 

Theorem 11 The derivatives of A{uj'^ , k'^) at anyuj'^,k'^ G C (hence in both the passbands 
and the stopbands at uj'^,k'^ G Mj are given by the formula 

WrkW^ = 2(-'' ' »""' /„ ^« X i (29) 

X F (<^i, . . . , M2 (?n+m + 1, ?l) M2 i<^l,<i2) ■ ■ ■ M2 , , 

where M2 {yi,yj) is a right off-diagonal component of the matricant M {yi,yj) , and 

F {<,!,..., ?n+m) = ^ /ai (?l) . . . fa„+m (Wm) , /o (?) = P (?) , fl (?) = 1^2 (?) ] 

CI = |(o-i, . . .,an+m) ■■ o-j = 0, 1; Q-t = m | , 

i.e. is a set of C^+m = (« + m)\/n\m\ permutations of a set (ui, . . . , cr„_|.m) , in which 
each (Ji is either or 1 and their sum is m. 

Proof. The expression ( |29| ) follows from the following property of matricants of related 
systems M- let Q(y)M(y,yo) = ^M(y,yo) and Q(y)M(y,yo) = ^M{y,yo) where 

Q(?/) = Q(y) + Qi(?/); then 

M (y, yo) = M (y, yo) /" [I + M (yo, ?) Qi (?) M (?, yo) d?] 

= M(y,yo)+ rM(y,?i)Qi(?i)M(?i,yo)d?i + ... (31) 

+ / d?i . . . / ' ' d<?jM (y, ?i) Qi (ft) M (ft , ?2) Qi (?2) . . . M (?j-, yo) + . . . . 
yo yo 



(30) 



Next note that Q [y; uj'^ , k'^) = Q [a;^,fc^] defined by (10)2 is linear in both cj^ and k'^. 
Denote small perturbations of and k'^ by e^j and e^. From (10)2, 



Q[oo^ + e^,k^ + ek]=Q[oo^k^]+i{f,2ek-pe^)T, T = (^J . 



(32) 
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Equation (31 ) with Qi = i {ii2£k — P^ui) T is therefore a Taylor series of M = M [w^ + e^^, /c^ + 
about the point £uj = 0, £k = 0, and hence the derivatives of the monodromy matrix M(l, 0) 
with respect to and /c^ are 

"1 r + m — 1 



a"+'^M(l,0) 



(-if i™n!m! 
X . . . 



d?i . . . 



(33) 



,) M (1, ft) TM (ft, ?2) r . . . M i^n+m, 0) 



with F defined in (30K Note that F = p (ft) . . . p ((j„) at m = and F = p2 (ft) • • • /U2 (?m) 

^trM(l,0) together imply 



at n = 0. Equation (33) and the definition A(a;^, k'^ 

gn+m^(^2^^2) 1 5'^+"^trM(l, 0) ( 



a(a;2)«a(A;2)™ 2 a(a;2)"a(A:2)' 



dft . . . 







+ m — 1 



+m 







(34) 



X F(ft, . . . , wm) tr [M(?„+„ + 1, ft)rM(ft, ?2)r . . . M(^„+^_i, ?„+™)r] , 



where we have used the identity tr [M (1, ft) . . . M (?„+m, 0)] = tr [M (?n+m, 0) M (1, ft) . 
and the fact that M {(^n+m, 0) = M (?n+m + 1,1) due to periodicity. By definition of T, 



Mr 



M2 
M4 



tr 



]V[«r...M('=)r 



.M. 



(k) 



(35) 



which reduces (34) to the desired form (29). 



Corollary 12 The first- order derivatives of A{uj ,k ) follow from (29) as 



dA 



p{y)m2{y)dy, 



dA 

djk^ 



P2{y)m2iy)dy, 



(36) 



where im2{y) = M2 {y + 1, y), see (15). 



Interestingly, the expression (29) for any derivative of A(a;^,/c^) involves, apart from 
piy) and/or P2{y), only a single, right off-diagonal, element M2(ft,?j) of the matricant. 



Recall that ReM2 = by (|13[)2, which conforms that (29) is real as it must be. Next we 
will obtain a different representation for the first derivatives of A(a;2, k'^) that is expressed 
via an eigenfunction u{y) of ([s]). In contrast to (29), this representation is restricted to 
the passbands |A(w2,A;2^j < i ^j^ji hence to co"^, k'^ G M. We note that the components of 
eigenvectors of M (1, 0), which appear in the explicit formulas below, are understood to be 
referred to a basis observing the identity (13) (an obvious counterexample is the Jordan 
form of M(1,0)). 

Theorem 13 The first derivatives of A{uj'^,k'^) within the open passband intervals A € 
(— 1, 1) (and hence uj'^, A;2 G satisfy the formulas 



dA 



sin K 
w+Tw 



piy) \uiy)fdy, 



dA 

djk^ 



sva.K 
'w+Tw 



P2{y)\u{y)\^ dy, (37) 



16 



where w is an eigenvector o/M(l,0) corresponding to the eigenvalue 



iK 



and u{y) is 



the first component of the vector ri{y) = M(y, 0)w = (u, ifiiu' 



At the cutoffs A = ±1, 



Eq. (31) yields zero derivatives in the exceptional case of a ZWS, and is otherwise modified 
to 



dA 



1 



2iwJTwg JO 



P{y) \u{y)\^dy, 



dA 



1 



2iwJTwg JO 



fi2iy)\uiy)\^dy, (38) 



where and Wg are the proper and generalized eigenvectors of M(1,0) that realize its 
Jordan form (see (44 V' ^'^'^ u{y) is equal to the first component of the vector rj{y) = 
M(y,0)wd. 



Proof of (31). The monodromy matrix M(l, 0) at | A| 7^ 1 has distinct eigenvalues q ^ q ^ 
and hence hnear independent eigenvectors wi, W2. Specify their numbering as 



M(l,0)wi = gwi, M(l,0)vi^2 = ^^^2 with g = e'^ / g 



-1 



iK 



(39) 



According to (31) and (32), 



- /'M(l.,)§<|lMfe,0)d!,= -iM(l,0) /'pWd,, 

Jo ) Jo 

^trP(y) =p(y)trr = 0). 



9(w2) 

where P{y) = p{y)M'^ {y , 0)TM{y , 0) 



(40) 



Hence, the derivative of A = ^trM (1, 0) at | A| 7^ 1 is 



dA 



1 

2i 



q [ Pu{y)dy+- [ P22{y)dy] =smKT [ Pii{y)dy 
.Jo Q Jo J Jo 



(41) 



where Pu is the upper diagonal element of P(y) in the base of vectors wi and W2. For 
the passband case A G (—1,1) being considered, the identity M^-*^ = TM+T (see (13)i) 
implies that 



wfTw2 = 0; Wj^Twi, w^Tw2 / [(w^Twi)(w^Tw2) < O] . 



(42) 



Using (|42j), the equality TM 
of r given in ( 32 ) , we find that 



-1 



(Mwi) T (following from (13)i) and the definition 



Piiiy) 



w+TP(y)wi _ p(y)77+(y)Tr?7i(y) _ piy)\u{y)\' 



Ae(-i,i) 



w^i^Twi 



w]^Twi 



w^Twi 



(43) 



where r]i{y) = M(y,0)wi = {u,ifiiu')'^ . Based on the numbering in (39) it follows that 
rji (1) = e*^wi and so u is an eigenfunction of (|5|) (see Corollary K3| . Substituting (43) 



into (41) and setting wi defined in (39) as wi = vi^ leads to (37)i. The proof of (37)2 
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is the same. Note that the sign alternation (23) of both derivatives at successive cutoffs 
is described in (37) by the factor (w+Tw) sini^ as follows: using K G [0,7r] implies 
s\nK > and alternating sign of w"*"Tw (due to switching between right- and leftward 
modes at successive cutoffs); while using unrestricted K > Q implies w^Tw < (rightward 
mode) and alternating sign of sini^. ■ 

Proof of (38). Consider a cutoff A = ±1 that is not a ZWS and hence implies a non- 
semisimple M(1,0). Denote 

M(l,0)wrf = gdWd, M(l,0)wg = grfW3 + Wd at A = = ±1, (44) 

which defines (not uniquely) the pair and Wg as a basis in which M(1,0) at A = ±1 
has upper Jordan form. Hence 



^A 

5(^ 



1 5M(1,0) 



1 

2i 



P2i{y)dy, 



(45) 



where P21 is the left off-diagonal of P(?/) at A = ±1 in the vector basis of and w^. The 
identity M"^ = TM+T for a non-semisimple M(1,0) implies that 



w+Twd = 0; w+Twg / [Re w+Twg = for det M = l] 



By (l46j) and the definition ^2 of P(y) 
P2i{y) 



w+TP(y)wd _ p{y)4{y)Trr^M _ p{y) Hy)Y 



A=±l 



wtTwfi 



(46) 



(47) 



where rjd{y) = M(y,0)wd = {u,iniu') . Inserting (47) in (45) provides (38)i. The proof 
of (38)2 is the same. ■ 



Note that (38) can also be obtained directly from (37) by taking its limit as |A| < 1 
tends to |A| = ±1. To do so, proceed from (39) with q, tending to Qd- It is al- 
ways possible to choose wi, W2 so that they have as a common limit and then 
(wi — ■W2) I {q — q~^) tends to w^, where w^^ and Wg satisfy (44). By using this limit- 
ing definition of Wg and the property w^Tw2 = (see (42 )i), the limit of the pre-integral 
factor in (37) with w = wi corresponding to q = e^^ is found to be 

sini^ " 



q-q ^ _^ 1 

2iwj^T(wi - W2) A^±i 2iwjTwg' 



(48) 



The factor wJ'Tw^ may also be expressed in terms of the elements Mj(l,0) = Mi of 
the matrix M(1,0) which satisfies (13). Using (44) yields two alternative forms of this 
expression as follows: 



w+Twg 



\Mi-qd\^ + \M2\^ 



IM4 



qdf + m\^' 



(49) 
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2qd, det [M(l, 0) - g^I] = and If Mi = qa hence M4 = Qd 

then either M2 = or M3 = 0, as occurs for instance if Q(y) is even 



If Ml, M4 7^ Qd then M2, M3 7^ 0, and so both formulas in (49) are equivalent, which 
follows from trM(l, 0) 
(or vice versa) 



about the midpoint of the period [0, 1], see the end of §3.1. Simultaneous vanishing of both 
M2, M3 is ruled out for a non-semisimple M(1,0). 



In conclusion, the combination of results (36) and (37), (38) yields the following inter- 
esting observation. 



Corollary 14 The right-hand sides of (36) are equal to those of (31) in the passbands 



A S (—1, 1), and to those of (38) at the cutoffs A = ±1 (unless the cutoff is a ZWS). 



3.5 Properties of the function m2{y) 



An important role of the function m2(y) defined in (15) is revealed by the fact that 



according to (36), the first derivative of A{uj'^, k"^) in or k"^ is an integral of m2{y) with 



a positive weight factor p{y) or ^2(2/)- Recall also that zeros of m2{y) are the Dirichlet 
solutions for the interval [y, y + 1], see §3.1. 

Theorem 15 The continuous function m2{y) = m2(y+l) satisfies the following properties: 
(i) if A{uj'^, k"^) G (-1, 1) then m2{y) has no zeros for y € [0, 1]; (ii) if A{uj'^, k"^) = ±1 then 
^2{y) > for any y E [0,1] or m2{y) < for any y E [0,1]; (Hi) if A{uj'^,k'^) ^ (—1,1) 
and Lip' , A;^ E M, then m2{y) has only finite number of zeros in [0, 1]. 

Proof. Consider (i). Suppose that A E (—1, 1) and there exists y such that 7712 (y) = 0. 
Then M (y + 1, y) has eigenvalues nii (y) and 771,4 (y)(= mrii^ (y) by det M = 1). Therefore, 
with reference to Remark [l| A = ^ \jni (y) + rn^^ (y)] , where mi according to (15) is real 



(since w , A; E M by Lemma ([5j)). Hence | A| > 1, which contradicts the initial assumption. 
The statement (ii) follows from (i) and the analyticity of A(ci;^,A;^). Consider (Hi). First 
note an identity 



M'(y+l,y) = Q(y)M(y + l,y)-M(y + l,y)Q(y) ^ 7T7'2(y) 



mi(y) - m4(y) 
Aii(y) 



(50) 



where ' = d/dy (if y is a point discontinuity of a piecewise continuous Q(y), then d/dy is a 
right or left derivative). Since //i(y) > 0, it follows that m^iy) = iff m/i{y) = mi(y). Now 
let us suppose the inverse of (Hi), i.e., that A ^ (—1, 1) admits the existence of an infinite set 
{yn}i° for which 7712 (y^) = 0. Without loss of generality we may assume that lim„_j.oo y-n = 
yo E [0,1]. Then m2{yo) = and 771,2 (yo) — 0- shown above, 777,2(yo) = yields 
7774 (yo) = (yo) and so we have A ^ (-1, 1) for A = | [777,i(yo) + ?77j"^(yo)] ^ (-1, 1) 



where ttii is real due to oj , k £ M. It therefore follows that 7774 (yo) 
According to (50), this contradicts 777,2 (^o) = 0. ■ 



mi ^ (yo) / mi (yo) 



The above result together with Eq. (36) provides a simple criterion for a ZWS, which 
complements Proposition [9| 
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Proposition 16 The following statements are equivalent: (i) (w, A:) is a ZWS; (ii) m2{y) = 
for any y. 



Proof. Assume (i). Then M(1,0) = ±1 by Proposition |9j Hence by (16) M(y + l,y) = ±1 
and so m2 = 0, which is (ii). Now assume (ii). It requires that A = ±1 by Theorem 
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dA{uj\k^)/d{ 



and yields dA/d{uj-^) = by Eq. (36)i. According to Proposition^ A{u}^,k^) = ±1 



impHes that {oJ,k) is a ZWS, which is (i). 



Interestingly, the function m3{y), whose zeros are the Neumann solutions for the in- 
terval [y,y + 1], shares some, but not all, of the properties of m2{y)- For instance, m3{y) 
displays the same properties (i), (ii) stated by Theorem 15 for m2{y) but it does not 
have the property (Hi). The dissimilarity stems from the fact that (50)i yields ra'^{y) = 
(a*2^^ — puP') (mi — 7714) , where, in contrast to (50)i, the first factor is not sign-definite. 
Also the derivatives of A(a;^,A:^) are not expressible via ms{y) as they are via m2{y) in 



(36). As a result. Proposition 16 does not hold for m3{y) in the sense that while it is true 
that m^{y) = for any y if {uj,k) is a ZWS, the inverse statement is not. An immediate 
counter-example is the point a; = 0, k = 0, where m'i{y) = for any y by (20) but this 
point is not a ZWS; moreover, the model case fJ-2{y)/p{y) = const = mentioned in §3.3 
ensures = on the whole cutoff line uji (0, k) = ck (see (|28[)) which has no ZWS points. 



Thus, the Dirichlet solution a;D,n(^) for [y, y + 1] does not depend on y only if {ujY),n, is 
a zero-width stopband, but the same is not generally true for the Neumann solutions. 



4 The dispersion surface /c) 

In this Section, we address the multisheet surface ojn {K, k) = ^u:\(K^ k) (> 0) which is 
defined by Eq. (18), and study the curves in its cuts taken at constant K, constant k and 
constant uj. 



Remark 17 If Eq. (18) with either K ork oruj being fixed defines a differentiahle function, 
then its derivative of any order can be expressed in terms of partial derivatives o/A(a;^, k"^) 



given in (29). 



Below we examine in detail the first non-zero derivatives. The higher-order ones are easy to 
obtain in a similar fashion by differentiating ( |18| ). It is understood hereafter that cj, A; G M. 
By (18), u}n{K,k) = ^n{—K-,k) = (K, —k) which permits confining considerations to 
ReK >0, k>0. 



4.1 The function UJn{k) for fixed K 



Consider the dependence oiconik) = uJn{K, k) for fixed K, Fig.[2j By Eq. (18), the branches 
ujn{k) are defined as level curves A(a;^, fc^) (= cosK) = const, which lie in the passbands for 
fixed E M 44> I A| < 1 and in the stopbands for fixed complex K G vrZ+i (M\0) \A\ > 1 
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(note that the branch numbering (22) does not apply in the stopbands, see the discussion 
of Fig. [2] below). 




Figure 2: (a) (left) The curves a;„(fc) = ujn{K, k) at different fixed K for the profile 
(21 ). (b) Sections of the curves for the piecewise constant profile defined by /ii = /_t2 = 
1, p = 1 for 2/ e [0, 1/2) and ^1 = ^2 = 12, p = 2 for y e (1/2, 1]. 



Proposition 18 If u / and d/S./d{uP') / 0, then 



k Aul _ k dA/d 
~dk ~ ^d(A;2) ~ ~^dA/d{uj^) 



(51) 
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where by and ^ 

_ /oV2(y)"i2(y)dy 



In addition, 



/oV2(y) \uniy)\'^dy 



K(^R or 



Iq p{y) \un{y)\'^dy 



(52) 



Km 



duji 
~dk 



dw„ 



LJ=0 

A:=0 



d/c 



k=0 



0, 



dyfc 
dwi 



ui=0 



(53) 



in view of 
7^ (note that uj 



The two other equahties in (53) follow from (51) and 
0, A; 7^ belongs to the stopband area where (52 )i applies, 



The former equality follows from (19) or else from (52) where m2{y) and ui{y) are constant 
at w, k = 

see Fig.[2}i). 

For G M, the excluded case dA/d{uj'^) = in (51 ) is related to ZWS discussed in §3.3. 
According to Proposition [oj if dA/d{uj'^) at X G M becomes zero then so does dA/d (A;^) 
and their simultaneous vanishing implies a ZWS. Barring extraordinary cases mentioned 
in 3.3.3, ZWS is an intersection point {oj, k)^^^ of two analytic curves Un (k) (as rigorously 
confirmed in Proposition [19] below), so there exist two derivatives at {uj, k)^^^. Their values 
can be determined by continuity from either of equations (52) applied in the vicinity of 
{uj, k)^^^ . Note that Eq. ( [52|i is not defined strictly at {uj, k)^^^ (where m2{y) = Vy, see 



Proposition 16) while Eq. (52)2 is, provided that Un{y) implies two different eigenfunctions 



from a subspace corresponding to two intersecting curves cj„ {k) at (w, k)^^^. 
Proposition 19 The curves ujn{k) for fixed K £M are monotonically increasing at k > 0. 



Proof. The function a;^(fe^) is analytic for any K G 
operators^of Kato's type A ^2\■ Hence if duj'^/d{k'^) 



since AK{k) is a family of analytic 
for some real k^, then there exists 



complex k^ in the vicinity of k'^ for which uj"^ = a;^(A;^) is real. But this would mean that 
the operator Bk{^) has a complex eigenvalue k'^ equal to k"^, which is impossible. Thus 
ujn{k) at ii' G M is a monotonic function. It increases by virtue of (52)2. To provide a fully 



self-consistent proof within the operator approach, note that (52)2 can also be obtained by 



applying the perturbation theory [15] to Ak given by ([6]), so that 



duJZ 



d(A;2) d(A;2) 



d {AKUn,Un), 



dAx 

d(F) 



Jq IJ.2{y) \un{y)\'^dy 
So Piy) l^n(y)|^dy 



(54) 



Consider the example plotted in Fig. [2| It demonstrates monotonicity of the curves 
ujn{k) = ujn{K,k) at fixed G M by tracing the cutoff curves at K = 0,11 {<^ \A\ = 1) 
and the curves at K = 7r/2 {-^ A = 0) within the passbands. Figure [2] also shows 
that, by contrast, the curves uj {k) = uj{K,k) in the stopbands, i.e. at fixed complex 
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K G vrZ + i (M\0) the level curves |A| = const > 1), may be not monotonic and can 
take a looped shape, either semi-closed or even fully closed. Note that the numbering 



of such curves cannot be defined by the rule (22) restricted to the passbands. A looped 
shape is due to a vertical tangent at a point where dA/d{uj'^) = (see (51), (52 )i). In 
any stopband except the lowest one, there exists a pair of curves a;cxt(^) and k^xti^), on 
which |A(u;^, A;^)| = cosh (ImK) has maxima in and in k"^ (in k aX k ^ 0), respectively. 
Hence each stopband except the lowest must contain looped curves uj (k) with a vertical 
tangent as they cross the curve LVext{k) - unless the latter fully merges with A;ext('^) as in 
the model case ^2{y)/ p{y) = const mentioned in §3.3. The curves c<Jext(^) and fcgxt (w) may 
intersect within a given stopband thus indicating a saddle point or an absolute extremum of 
A((x'^, k"^) (the latter is exemplified in Fig.jSja, see the family of closed level curves |A| > 1). 
At the same time, oj^^tik) and A;ext('^) cannot contact the cutoff curves except at the point 
of a ZWS (see Fig. ^)), which is always a saddle point of A((jj^, k"^). 

It is shown in Appendix A3 that the lower bound for the branches a;„(/c) at IT G M 
is miuj^gjo,!] a/zWP- In the remainder of this subsection we prove that this bound is also 
a common limit of ojn{k). To do so, it is convenient to introduce the velocity u„ = oJn/k. 
First we specify the derivative of Vn{k) in order to demonstrate its monotonicity (note that 
it is easy to similarly obtain sign-definite derivatives at fixed X € M for any other optional 
choice of the pair of spectral parameters among w, k and v oi s = v^^). 

Lemma 20 Let K ^ n £ N be fixed. Then f^(A;^) = (fc^) /k'^ is a decreasing function 
with derivative 

dvl ^ 1 /oVi|<(j/)pdj/ 



where Un and u'^ are defined by ri{y) = (n, ifiiu') - = M(y,0)w taken at col (cf. (S^) 



Proof. Multiply Eq. ([2]) by u (= Un), integrate by parts and divide the result by fc^, to 
yield 

j p\un{y)\'^dy = ^ ^i|n„pdy + j p2\un\'^dy. (56) 



Substituting from (56) along with (54) into da;^/d(A;^) = k'^dv^/d{k'^) + leads to (55 
The same result follows by applying the perturbation theory |15) similarly as in (54) 
whence dT;^/d(A;^) = —{{piu'^)',Un)p/k'^ ll'^^nllp and integrating by parts yields (55). ■ 

Proposition 21 Let K be fixed. Then for any n gN 

hm 4 = min (57) 
A:-s>oo s/e[o,i] p{y) 
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Proof. Rewrite ^ in the form 



/f2 

p 



0. 



(58) 



where v'^ = uP' jlP'. For any fixed v = a> min ^J 112/ p, the coefficient (/i2/p) — v'^ changes 
sign on the interval [0, 1] and hence there exist infinitely many distinct values fc^ > which 



satisfy (58) (see more in [9J). The latter means that any curve Vn{k), n E N, intersects the 
line a{k) = a for any a > min \J p. Combining this statement with the above-mentioned 
facts that all Vn{}i) are decreasing and have the lower bound min \J 112I P yields (57). ■ 



It is noteworthy that there is no common limit for a finite spectrum of eigenvalues of a 
discrete Schrodinger operator with a large potential [E 
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Figure 3: The Floquet branches LOn{K) = uJn{K,k) at fixed fc = 1. (b) Real isofre- 
quency branches Kj{k) at fixed oj ~ 8. The same profile (21) is used. The cutoff 
values of uj in (a) and of k in (b) can be compared with Figs. |l] and [2^. 



4.2 Function uJn{K) for fixed k 



Consider the function ujn{K) = {K, k) implicitly defined by Eq. (18): A(tj^, k'^) = cosK 
at fixed k. Since iOn{k) is periodic and even, it suffices to deal with one-half of the Brillouin 
zone Reii' G [0, vr], see Fig. [3^. For brevity, denote the cutoff values u;„ (vrm, k) of ujn {K, k) 
as 

UJniTTm,k) = UJn,m, m = 0,l. (59) 



Let us indicate the passbands and stopbands of u}n{K,k) by ImK = and ImK ^ 0, 
respectively (the latter being short for K = vrm + i Im K ^ vrm) . Explicit expressions for 



the first non-zero derivative of uJn{K) readily follow by expanding both sides of (18) and 
invoking the formulas for dA/d{co'^) obtained in §3.4. Note that Eq. (60) with (37)i for 
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real K (see below) can also be obtained by means of perturbation theory |15j applied to 
an appropriately modified form of ([2]), ([3]) with an operator explicitly dependent upon K. 

Proposition 22 If either (i) Imi^ = and K 7^ vrm (hence dA/uj ^ by Proposition^ 
or (a) Im K and dA/duj / 0, then 



sinK 



dK (dA/duj) 



(60) 



where sin K = Vl - A2 and dA/du = 2ujdA/ diu"^) is given by aSm) or ([5?J)i for (i) and 

m an 

dUn 



by (36)i for (ii). If K = irm and dA/doo 7^ 0, then 

(-1)™+^ 



0, 



d^uj„ 



(61) 



dK dK^ idA/dio),, ' 

where ujn,m = ^n,m{k) are the roots of equation A(u;^, /c^) = (—1)™ and dA/duj is given by 
or'^i. 

Consider the special cases where dA/doj = 0. Let K = vrm and dA/du = at cj / 0, 
which implies a cutoff ujn,m corresponding to a ZWS. Then 



dcJn 

dK 



(-l)-+"+7^(-l)-+i(a2A/5u;2)^_. 



(62) 



Next let ImK 7^ and dA/doj = 0, which defines the point uj = Woxt in a stopband at 
which |A(a;)| = cosh(ImX) reaches its maximum |Aext| > 1 (see Fig. [2] and its discussion 
in §4.1). The function Imi^(a;) satisfies [dim. K / duj) = and 



d^ IraK 

duj"^ 



(a^A/a 



a/a 



2 

cxt 



1 



;< forlm/C > 0) . 



(63) 



The explicit form of d'^A/duP', which appears in (62), (63) and is negative at m = and 
defir 

92 A 



positive at m = 1, is defined by (29). It can be written in the following equivalent forms 

d'^A 



a(L<j2)2 
1 



dy 







-4w 

y+l 



dy / p{y)p{yi) M2{yi + l,y)M2{y,yi)dyi 
Jo 

p{y)p (yi) M2 {y + 1, yi) M2{yi,y)dyi (64) 



y 



-2iJ I dy piy)piy + yi)M2{y + l,y + yi)M2{y + yi,y)dyi, 
Jo Jo 
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where dA/dco = and u) ^ (i.e. dA/d{uj^) = 0) have been used. Fmahy, consider the 
case CO = 0, which imphes dA/duj = 0, d'^A/duJ^ = 2dA/d{uj'^). If both a; = and k = 



K = 0), then referring to (19), the derivative (62) for m = 1 reduces to 



^ = ^/^ipUjH^y (65) 



If a; = and A; > K = ilmi^ / 0), then (dlmK/da;)^^^ = and (63) becomes 



d^Imi^ _ 2[dA/d{oj 



2\ 



Jw=0 



da;2 {0,k^) - 1 



(66) 



where [dA/d{uj'^)]^^^ < is given by (36) 



It is evident from Eq. (60) that the Floquet branches u)n{K) for any fixed real k are 
monotonic in G [0, vr]. For completeness, let us also mention two important results from 
the general theory of Schrodinger equation [T^l [ISl E] that extend to the case of Eq. ^ 
with fixed k. These results state that Imi^(a;) is a convex function and that each branch 
u}n{K) has one and only one inflection point in K G [0,'7r], unless it is the lowest branch 
uJi{K) at A; = or a branch bounded by a ZWS at either or both cutoffs K = vrm, in which 



case there is no inflection points. Note in conclusion that Eqs. (|61j) and (62) provide an 
explicit definition for the near-cutoff asymptotics of branches ujniK) that were analyzed in 
[7] by a different means (the scaling approach, also extended in [7j to 2D-periodicity). 

4.3 The function K(k) for fixed u 

Consider the dependence of K[k) = arccos A{uo'^, k"^) on /c > at fixed cj. Let the branches 
Kj[k) £ [0, vr] for real K be numbered in the order of increasing k. Since w„ (k) = uJn{K, k) 
is strictly increasing in k (see Fig. [2]), the number of real branches Kj{k) at any fixed value oj 
is fully defined by its position with respect to the frequency-cutoff points at k = 0: there is 
a single real branch Ki{k) for a fixed co in the interval < < ^2 (vr, 0) ; two real branches 
Ki{k), K2{k) for u in u}2 (tt, 0) < oj < a;3 (0, 0) ... etc. Besides, the first real branch Ki (fc) 
starts at A; = and spans a range [0, vr) or (0,vr] iff |A(a;2,0)| < 1, i.e. iff the given u is 
fixed within the passband at A: = 0. For example, the value w = 8 G (a;3(0, 0), Li;4(vr, 0)) in 
Fig. [2] yields three real branches Kj{k) with Ki{k) G [0,vr), see Fig. [sja. 
Denote by 

kj,m{uj) = kj^rn, m = 0,1, (67) 

the roots of equation A{uj'^, k"^) = (—I)™' which define the points at which Kj{k) = vrm and 
the given uj is the cutofFj these points kj^^n 

are separated by the stopband intervals |A| > 1 
where ImK ^ 0. The explicit form of the first derivative of K (k) for real or complex K 



follows from (18) and the formulas for dA/d{k'^) in exactly the same way as that ujn{k) in 
U.2. 
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Proposition 23 If K ^ vrm and k 0, then 

dK dA/dk 



dk 



sinK ' 



(68) 



where dA/dk ^ for real K. If Kj (k) = vrm at /c / and {dA/dk)j^, ^ ^ 0, then the 
locally defined inverse function k (K) satisfies 



dk d'^k 
0, 



\m+l 



dK 



dK'j {dA/dk)^ 



(69) 



If k = 0, then 



dK^ _ \ 0, ^ = -sii^ {dAld{k ^)\ at K^ / vrm, 

at K\ = vrm. 



dA: 



2(_i)"+i [aA/9(fc2)] 



(70) 



fc=0 



Consider the imphcation of possibly existing ZWS. Assume that a fixed u is a ZWS for 



some k ^ 0. This means that Kj (kj^m) = vrm and {dA/dk)^ 



where kj^m 7^ 0. Then 



( 69 ) is altered to 



d^ 
dk 



(-1) 



m+j 



(71) 



Now assume that a fixed a; is a ZWS at A; = 0, i.e. let Ki 



■Km 



Then dKi/dk = by (70), and 



and [dA/d{k'^ 



]fc=0 



d^Ki 
dk'^ 



2(_i)'"+i [S2A/9(fc2) 



Jfc=o- 



(72) 



The second-order derivative of A in (71), ( |72| ) can be obtained by differentiating (86)2 in 
the same way as in (lei. Note that d^/d{9f also appears in the formula analogous to 
(63) for d2lmi^/dA;2 at the point k^xt where dlm.K/dk = 0. 



Thus, by (69) and (71), all real branches Kj{k) at fixed cj have vertical tangents at the 
edge points Kj {kj^m) = ^rm, kj^m 7^ (see Fig. [ijs), unless the cutoff io = ujn (7rm, kj^m) is 
a ZWS in which case Kj{k) does not make a right angle with the line K = vrm. In turn, 
by (70) and (72), the real branch Ki{k) has a horizontal tangent at k = 0, K ^ tt and a 
non-zero first derivative at A; = 0, X = vr, unless u = iVn (t^, 0) is a ZWS stopband in which 
case the slope of Ki{k) vanishes at A; = 0, X = vr. 

Remark 24 If the cutoff uj = oon (vr, 0) is not a ZWS, then (i) the curve Ki(k) = Ki (—A) 



has a kink at k = 0; (ii) \/uj{K, k) = at k = by virtue of (53) and (61) 
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4.4 Convexity of the closed isofrequency branch Ki (k) 

The normal to real isofrequency branches Kj{k) defines the direction of group velocity 
Vuj{K^ k) which makes their shape relevant to many physical applications. In particular, 
negative curvature of an isofrequency curve is known to give rise to rich physical phenomena 
related to wave-energy focussing. Since the function K{k) = arccos A with |A| < 1 defines 
a unique K G [0, vr] , no vertical line can cross twice the curve K(k); however, this by itself 
does certainly not preclude a negative curvature. In fact any real branch Kj{k), which 
extends from Kj = to Kj = vr, has vertical tangents at those edge points and hence must 
have at least one inflection between them (unless the exceptional case of ZWS, see §4.3). 
This simple argument, however, does not apply to the first branch Ki{k) if the reference 
UJ is taken within the passband range at fc = and hence Ki{k) does not reach one of the 
edge points or vr. In other words, the situation in question is when Ki{k) extended by 
symmetry to any real k^Q forms a closed curve. 

In the present subsection we address an important case of a relatively low frequency 
w which is restricted to the passband below the first cutoff wi (vr, 0) at the edge of the 
Brillouin zone X = vr at /c = 0. For any fixed uj < (vr, 0), there is a single real isofrequency 
branch Ki{k) = arccos A (cj^, /c^) G [0,vr) that is continuous in the definition domain k G 



[— ^i,0;^i,o] ) where ki^Q is the least root of equation A = 1 (see (67)). According to (91 )i, 

(p) /(M2) < kifl{uj) < wmaXj^g[o,i] \/p{y)/My)- (73) 
We will show that Ki{k) is strictly convex. The proof is preceded by a lemma. 

Lemma 25 For fixed u < oJi (vr, 0) , derivatives of the function A(u;^, k"^) of any order in 
k"^ are strictly positive at k"^ > 0. 



Proof. Let uj = {). Then A(0,A:^) > for A;^ > by (l84j) and so dA{^,k^) / d{k^) > for 
because A(A;^) at fixed o;^ satisfies the conditions of the Laguerre theorem (see 
Proposition [7| . In other words, all zeros of dA{0, k'^)/d{k'^) lie in k^ < (see Fig. ^p)- 
Now let < a; < a;i (vr, 0) . This means that —1 < A(w^,0) < 1 and so the first zero of 
5A(a;^, fc^)/c?(/c^), which is where A < —1, still lies in k"^ < 0. Thus, if u; < wi (vr, 0) then 
5A(cj2, A;2)/a(/c2) > for > and hence, again by the Laguerre theorem, dPA/d{k'^)P > 
for A;^ > and for any p > 1. M 

Theorem 26 The curve Ki (k) is convex at any fixed oj such that uj < uji (vr, 0) . 
Proof. The second derivative of Ki[k) is 

K'l{k) = -{1-Ar"'h, h{k)^AQ)\{l-A^)^^, (74) 



where -1 < A^ < 1 for A; € (-A;i,o, A;i,o), see (jeTj). Note that dA/dk = at /c = 
0. Let CO < uji{tt,0). Then h{0) = {l - A'^) d^A/dk"^ > and h'{k) = {dA/dkf + 
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(l - A^) d^A/dk^ > according to Lemma [25] Due to h{0) > and h' (k) > at k > 0, 
it follows that h{k) > at > 0. Hence K'^ik) < in its definition domain [— fci^O) ^i,o]- 
Thus, Ki{k) is convex. ■ 

The obtained result sets an important benchmark against any artefacts of approximate 
analytical and/or numerical modelling of the first isofrequency curve Ki{k) = arccosA, 
which are possible as a result of truncating series for arccos or for A = itrM(l,0) see 
(12)). Figure |4] demonstrates an example where an approximate computation of Ki{k) 
produces a spurious concavity. In this regard we note that Figure 1 of [20j, which is 
sketch of the generic relation between K and k for fixed but small ui, incorrectly gives the 
suggestion that concavities can occur. 
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Figure 4: (a) The approximate and (b) the exact first isofrequency curve Ki{k) = 
arccos [itrM(l, 0)] at fixed a; (= 3.4) < wi (tt, 0) for a periodically piecewise constant 
profile defined by /xi = 1, /12 = 0.35, p = 0.2 at y e [0, 1/2) and /xi = 0.95, 112 = 0.4, 



p = 0.19 at 2/ G (1/2, 1]. The monodromy matrix ( 12 ), which in this case is M (1, 0) = 
(exp Q2) (exp Qi) with Qj defined by ( 10 12, is computed via the series of the co-factor 



exponentials, keeping four terms for each of them in the case (a) and 30 terms in the 
case (b). 



In conclusion, a remark is in order concerning the high-frequency case where the first 
isofrequency branch Ki{k) defined in k £ [—^1,0)^1,0] is accompanied by the higher-order 
branches i^j>2(A;). In general, Ki{k) should stay convex and Kj>2ik) should have not 
more than a single inflection point. However, it seems possible to construct a theoretical 
example, though quite peculiar, of a periodic profile, for which the above is not true. 
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Appendix 

Al. Properties of the operators Ak and Bk 

It is evident that the operator Ak defined in ^ is symmetric for k'^, i^T £ M, i.e. 

{Aku,v) = — / (^fiiu'y v*dy + / fi2uv*dy = / niu' v'* dy + k^ / fi2uv*dy 
Jo Jo Jo Jo 

= -/ {fiiv')'* udy + k"^ fi2UV*dy = {u,AKv) , 
Jo Jo 

using the identities fiiu'v* |q= niv'u* |q= which follow from the boundary condition ([t]) 
on «, G Dk iff K is real. The proof of the symmetry of Bk for uP' , K the same. 

We now demonstrate that Ak and Bk are self-adjoint with discrete spectra a {Ak) = 
{o;^}^ and a{BK) = {^n}]" corresponding to complete sets of eigenfunctions (as stated 
in §2). This is achieved by explicit construction of the resolvent of each operator, TZk^x = 

(^Ak — oj'^) ^ or TZk,\ = {Bk — k'^) ^ , where A implies or k'^. In order to do so consider 
the equivalent equations 



(75) 



%lli]r4: itX:] ^^t.u{y)eDK,y{y)eLM, 



(76) 



which can be recast as 



v'{y) - Q{y)v{y) = i{y) with 7(y) = (.^^^y^) , (i) = e'^r?(o), (77) 

where / = —ipg for Ak, f = for Bk, and 77, Q are defined in ([T]), ([To]), respectively. 
The solution to (77) is a superposition of its partial solution rjp with the solution r/o(2/) of 
the corresponding homogeneous equation: 



v{y) = vp{y) + m{y), Vp{y)= M (y,?) 7 (?) d?, r?o(y) = M(y, 0)770(0). (78) 

Jo 

The vector 770(0) is found from the quasi-periodic boundary condition that yields 77^ (1) + 
770 (1) = e*^77o(0). Thus 



ri{y) 



[ G(77,?)7(?)d? with 

Jo (79) 
G {y, = M {y, ^) H {y - ^) - M{y, 0) [M(l, 0) - e^^l] M (1, ?) , 



where H {y — q) is the Heaviside function and e'"^ is not an eigenvalue of M(l, 0) for the 
given ^ a {Ak)-, k^ ^ o {Bk)- It can be checked that the Green-function tensor G (y, ?) 
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satisfies the identity G {y, <j) = 
satisfies Gi2 (y,?) = -Gl^ (y,?)- By 



TG^(?,y)T, so that its right off-diagonal component 



u 



T^K,\g = G{y, A) / (?) d?, where G (y, A) = iGi2 
Jo 



(80) 



It is seen that the resolvent TIk,x is an integral (bounded) self-adjoint operator generated 
by a piecewise continuous kernel. The symmetry {JlK,\g-,v) = {g,TlK,\'v) follows for any 
V £ Dk from G {y, A) = G* (?, y; A) or else from the symmetry of Ak, Bk- Thus TZk,\ 
satisfies the Hilbert-Schmidt theorem and Ak, Sk therefore possess the above-mentioned 
properties. 



A2. Bounds of the function A(c<;^,A;^) 

The far-reaching properties of the analytic function A(uj'^, k?) stated in Proposition [t] follow 
by applying Laguerre's theorem to A(a;^) at any fixed k"^ and to A(A;^) at any fixed w^. 
A function satisfying Laguerre's theorem must be an entire function of order of growth 
less than 2. Verification of this condition for A(a;^,A;^) requires its uniform estimation 
in C. The WKB asymptotic expansion (see §3.2) is not well-suited for the task in hand. 
Here we derive explicit bounds which show that A(a;^) and A(A;^) for, respectively, any fc^ 
and ijp are entire functions of order of growth ^. The derivation consists of two Lemmas 
in which the following auxiliary notation is used: /max = niax/(?/), /min = min/(y) for 
f{y) = p{y), l^i,2{y) and y G [0, 1]. 



Lemma 27 For any uj, k £ C, 



|A(a;^,A;^)| < cosh Y'/iimin(M2max + Pmax 

Proof. For any 2x2 matrix A with the entries (ai-.a^) , define |A| as 



Oil 
as I 



ml 

1 04 



(81) 



(82) 



and note that \Y\„ A„| < | A„| where the entrywise inequality is understood. Recall that 
/ appearing in (12) implies a product integral and is an exponential when the integrand 



matrix is constant. Hence it follows from (|10|)2, (|12|) and (|17|) that 
|A(^j2,A;2)| = 



tr/o [I + Q(y)dy] 



1 + i 



I + i 





max 



\UJ\ 



A^lmin 









f^2{y)k^ - p{y)uj^ 



-Pi\y) 




dy 



dy 



cosh 



max I 



Mil 



(83) 
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The inequality (83) confirms that A(w^) and A{k'^) are entire functions with order of 
growth not greater than ^ in each argument. Next we demonstrate that A for certain 
w^, k"^ grows no slower than an exponential of power ^ of and/or k"^. This will enable 
us to conclude that the order of growth of A(a;^) and A(/c^) is precisely ^. 

Lemma 28 For uj"^, k'^ e M, 



|A(a;^A;2)| > cosh y //^^^(^ammfe^ - Pmax^^^) for fc^ > /iaminPn 

Proof. First introduce a class Ai of 2x2 matrices such that 



M 



ai 
ids 



-ia2 
a4 



a,- > 0, j = 1..4. 



(84) 



(85) 



For two matrices A and B from A^, we say that A >_v( B iff aj > bj for any j = 1..4. If 
A e M and B e then AB e M also. Therefore, if A^, B^ G and A^ >m 
for any k = l..n then Ai..A„ '>m Bi..B„ and tr (Ai..A„) > tr (Bi..B„) (which is easy to 
check for n = 2 and is therefore valid for any n). We note from (10)2 that ^2min^^ > Pa 
implies I + Q(y)di/ G M. for any y G [0, 1] and dy > 0; moreover, 



l + Cl{y)Ay>M I + J 



^^2 min 



A'l max 





dy 



(86) 



and consequently 



A(t<;2,A:2) 




^ I + Q(y)dy] > ^tr 



.1 . 



max 





dy 



cosh 



- 

/ill 



(87) 



A3. Bounds of the first eigenvalue uf{K,k) 

Proposition 29 For K G [— tt, vr] and G M, i/ie /irsi eigenvalue ujf {K,k) is bounded as 
follows 

k mm — -— < LJi (A, /c) < ^^if + -r^k . (88) 
?/e[o,i] p{y) {p) {p) 

Proof. Let ui G -Div- with the unit norm = 1 be the eigenfunction vector of Ak 

corresponding to the eigenvalue w^. Then 

uj\ = (Akui.,ui) = / ^1 I n'l P dy + A;^ / /^2 I'^^il^ dy > A;^ / — p\ui\^ dy > k^ min — . 

^ 7o Jo Jo P J/6[o.i] P 
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An equivalent proof of the lower bound ( 89 ) follows by noting that the initial equation ^ 
yields zero as the sum of the positive operator — (/xiu')' and the operator multiplying u 



by 



implying that the latter factor must be negative. In order to obtain the 



upper bound, introduce the function v{y) = (p) e^^^ such that v{y) G Dk and 
Hence cof as a minimal eigenvalue of Ak satisfies 



u;^ 



inf 

ueDK, ||«||p=i 



{Akui,ui) < {Akv,v) 



ip) 



1. 



(90) 



Corollary 30 The hounds of the first cutoff at the centre and the edge of the Brillouin 
zone are, respectively, 



k min 

?/G[0,l] 



'My) 
p{y) 



< uji (0, A;) < k^ 



'{P2 



(P) 



uJi (0, k) < uji (vr, k) < 



(P) 



-7r2 + 



{P2] 
(P) 



■k^. (91) 



As stated in Proposition 21, the lower bound (88) of u}i{K, k) and hence of all curves 
ujn {K, k) for G M is also their limit at A; — )• oo. Note that oji (0, k) > uj-^^i(k) by (24), 



where a;N,i(A;) is the lowest branch of solutions of the Neumann problem for y € [0, 1]. It 
has the same bounds and the same limit at /c — t- oo as wi (0, k). In this regard, recall the 
model example P2{y)/p{y) = const = (see §3.2), where wi (0, A;) = a;N,i(A;) = ck merge 



together with their upper and lower bounds. By (pjl|)i, unless wi (0, A;) is a straight line, it 
has an inflection point (and so does a;N,i(A;)). Furthermore, the case of constant p, pi^2 is 



an elementary example of the equality of the upper bound in (88) and (91)2 
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